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ABSTRACT 

We derived the luminosity function (LF) of dwarf galaxies in the Coma Cluster down to Mr = — 10 
at three fields located at the center, intermediate, and outskirt. The LF (—19 < Mr < —10) shows no 
significant differences among the three fields. It shows a clear dip at Mr ~ —13, and is composed of 
two distinct components of different slopes; the bright component with — 19 < Mr < — 13 has a flatter 
slope than the faint component with —13 < Mr < — 10 which has a steep slope. The bright component 
(—19 < Mr < —13) consists of mostly red extended galaxies including few blue galaxies whose colors 
are typical of late-type galaxies. On the other hand, the faint component (—13 < Mr < —10) consists 
of largely PSF-like compact galaxies. We found that both these compact galaxies and some extended 
galaxies are present in the center while only compact galaxies are seen in the outskirt. In the faint 
component, the fraction of blue galaxies is larger in the outskirt than in the center. We suggest 
that the dwarf galaxies in the Coma Cluster, which make up the two components in the LF, are 
heterogeneous with some different origins. 

Subject headings: galaxies: clusters: individual: Coma Cluster (Abell 1656) — — galaxies: luminosity 
function galaxies: dwarf 



1. INTRODUCTION 

The galaxy luminosity function (LF) is a useful mea- 
sure to describe fundamental statistical properties of 
galaxy populations and serves as a clue to study the his- 
tory of galaxy formation an d evolution. Since the mon- 
umental work of the LF by Sandag e et al.l (|1985f ). rela- 
tively bright part of the LF ha s been studied extens ively 
both for cluster galaxie s (e.g.. iPopesso et al.l I2005T ) and 
for field galaxies (e.g., iBlanton et al.l 120031 ). However, 
our knowledge of the faint part of the LF (Mr > —19), 
which is dominated by dwarf galaxies, is still very lim- 
ited. This is due to the difficulty in determining the 
distance of a large number of intrinsically faint galaxies. 
Published results of the faint part of the LF, especially, 
those of the faintest part (Mr > —13) remain scarce. 

Clusters of galaxies are a good target to investigate 
the faint end of the LF because some statistical meth- 
ods can be applied to identify member galaxies. In this 
study we chose the Coma cluster in order to investigate 
the faint part of the LF, that is, to study properties of 
dwarf galaxies and their possible dependence on different 
environments in this cluster. The Coma Cluster (Abell 
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1656) at z = 0.023 (|Struble fc Rood|[l999l) is the richest 
cluster of galaxies in the local Universe. The virial radius 
(?"20o) of this cluster is 2.0 h~ x Mpc, which corresponds 
to a virial mass of 1.9 x lO ^h^Mp, as estima ted from a 
weak lensing measurement ()Kubo et al.l l2007) . It is rel- 
atively easy to take deep images of this nearby cluster 
which enables us to reach the faintest part of the LF. 

The LF of th e Coma Cluster was st u died by many 
authors (e. g.. IBernstein et al.l 119951: iBTviano et al 
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2002; 



Seeker et all 119971: iTrenthaml 119981: lAdami et al 
Andreon fc Cuillandrd 20021 : IBeiiersbergen et al 
' Mobasher et al.l [20031 iMilne et all 120071: 



Adami et al.l 120081 120091) . Several techniques were 
used to construct the Coma LF. The statistical back- 
ground/foregrou nd subtraction is historically the 
mam one (e.g.. IBernstein et al.l 119951: iTrenthairj [1998; 
[ Andreon fc Cuillandrd 120021: IBeiiersbergen et all l2002i 
IMilne et al.1120071 ). O thers include the spectroscopic red- 
shift technique (e.g.. lAdami et al.l 120001 : iMobasher et al.l 
I2003t lAdami et al.l|2009D photometric redshift technique 
(e.g., lAdami et al. | 20081) and color magnitude relatio n 
technique (e.g.. iBiviano et al.l 119951 : ISecker et al.l 119971) . 
The faint-end slope of the LF derived from previous 
studies do not agree with each other and it seems that 
it depends on the magnitude ranges where the slope 
was measured, i.e., the slope is flatter at brighter mag- 
nitudes, and it becomes steeper at fainter magnitudes. 
However, a simple comparison of the values of the 
slope may be inappropriate because of the difference in 
the observed area and galaxy sampling criteria among 
previous studies. Most previous studies observed the 
cluster core only, with a few probed outskirts as well as 
the core. 

The aim of the present study is to investigate the be- 
havior of the faint part, especially, the faintest part of 
the LF of Coma cluster using a statistically significant 
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large sample constructed from deep and wide photomet- 
ric images obtained with the Subaru Suprime-Cam. Our 
imaging data cover much wider area than most of the pre- 
vious studies that reached similar limiting magnitudes. 
We employ a careful procedure for subtracting the back- 
ground galaxy population statistically to derive the LFs 
of the cluster galaxy population. 

We investigate the properties of the very faint dwarf 
population and their possible environmental dependence, 
if any, in terms of LF, color, and surface brightness. Our 
galaxy catalogs for the core and the outer areas are con- 
structed by using the same instrument, same niters, and 
the same selection criteria, which ensures coherent and 
homogeneous photometric data to allow credible com- 
parisons of galaxy properties in different environments 
of the cluster. 

We adopt the standard cosmology model with Hq = 70 
km s _1 Mpc -1 , Q m = 0.3, and fl\ = 0.7. Magnitudes 
are on the AB system. We assume a distance to the 
Coma Cluster of 100 Mpc and a distance modulus of 
35.0 throughout this paper. 

2. OBSERVATIONS AND DATA 

2.1. Imaging Data 

We perform B- and J?-band observations with 
Suprime-Cam (Miyazaki et al. 2002) mounted on the 
Subaru Telescope (jive et al.l I2004f )for a total of three 
nights, between 12 and 14 May 2007. Observing condi- 
tions were mediocre due to the occasional cloud traver- 
sal on all three nights. The Suprime-Cam covers a 
34' x 27' field of view with a pixel scale of 0".202 pixel - 1 . 
A typical seeing size was ~ 0."8 in B- and i?-band. We 
selected three fields in the Coma Cluster for the present 
study, and they are hereafter referred to as the Coma 1, 
Coma 2, and Coma 3 fields, respectively. The Coma 1 
field is located in the central region of the Coma Cluster. 
Coma 2 is located at 1.3 Mpc projected distance to the 
southw est of the cluster core, in the postulated infalling 
region (jBiviano et a l. 1996; Neumann et al.|[2001[ ). Coma 
3 is about 1.8 Mpc away from the cluster core, i.e., the 
outskirts. Figure Q] shows the three observed fields on 
the digitized sky survey (DSS) image. The distribution 
of X-ray emission is shown in Figure [5J We summarize 
the details of observation of the Coma Cluster in Table 

m 



2.2. Control Field 



In order to obtain the intrinsic LF of cluster galaxies, 
we must correct for the effect of contaminant galaxies, 
which are mainly located behind the cluster. To esti- 
mate the number of contaminant foreground and back- 
ground galaxies, we use the Subaru Deep Field (SDF; 
IKashikawa et al.|l2004D 7 as the control field, which is lo- 
cated near the North Galactic Pole. The SDF was care- 
fully chosen to avoid bright stars, bright galaxies, and 
nearby clusters. This field was observed with Suprime- 
Cam and multi-band data based on the common seeing 
size of ~ 0."98 are available. We use the B- and i?-band 
data taken from the SDF archive. Their filter sets are 
identical to those we used for the Coma Cluster. SDF 
data reach much deeper limiting magnitudes than the 
Coma data. 

7 http://soaps.nao.ac.jp/sdf/Project. html 
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Fig. 1. — Digitized Sky Survey (DSS) (1st) image of the Coma 
Cluster. North is up, and East is to the left. The rectangles show 
the Coma 1 (center), Coma 2 (NGC 4839 groups), and Coma 3 
(outskirt) fields. The three giant elliptical galaxies are marked 
with circles, NGC 4889, NGC 4874 and NGC 4839. 
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Fig. 2. — ROSAT X-ray image of the intracluster gas in the 
Coma Cluster from NASA HEASARC. North is up, and East is to 
the left. The color scale shows higher X-ray intensity for red and 
lower one for blue. 

We consider that the population of Galactic stars in 
the SDF is similar to that in the Coma Cluster because 
the two fields have almost the same Galactic latitude 
(b ~ 88° for Coma, b ~ 83° for SDF). Accordingly, we 
use the SDF data to correct for the effect of Galactic stars 
in R > 22 where we do not apply star/galaxy separation. 

We also use the imaging data of the Subaru XMM- 
Newton Deep Survey fields (SXDS: lFurusawa et al.l2008j) 
8 to examine the effect of the cosmic variance (see £14.31 
for details). The SXDS fields are also observed with 
Suprime-Cam and the same photometric calibrations as 
the SDF was applied. 

3. DATA REDUCTION AND CATALOG 

3.1. Data Reduction 

Data reduction was carried out in a standard proce- 
dure usi ng the reduction software for Suprim e-Cam; SD- 
FRED (jYagi et al.ll2"002at lOuchi et al.ll200l and Image 
Reduction and Analysis Facility (IRAF) . We perform the 

8 http://www.naoj.org/Science/SubaruProject/SDS/datapolicyJ.html 
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Region 




Date 




R.A. 

(J2000.0) 


Dec. 
(J2000.0) 


Band 


Exp. Times 
(min) 


Lim. Mag. 
(5<t) 


Effective Area 
(aremm ) 


Coma 1 




12 Mar. 


2007 


12 h 59 m 46 s 


+27°54'02 " 


B 


60 


26.0 


769 




13 


14 Mar. 


2007 






R 


35 


25.9 




Coma 2 




12 Mar. 


2007 


12 h 57 m 22 s 


+27°20'02 " 


B 


60 


26.5 


880 




12 


13 Mar. 


2007 






R 


90 


25.9 




Coma 3 




12 Mar. 


2007 


12 h 57 m 23 s 


+26°53'09" 


B 


60 


26.3 


863 




12 


14 Mar. 


2007 






R 


114 


26.3 





TABLE 1 

List of the observations for the Coma Cluster. 



point spread function (PSF) equalization by convolving 
gaussians to the reduced Coma images in order to match 
their PSF size (FWHM) to that of the SDF image. The 
resulting PSF size is 4.9 pixels (0."98) in all the images, 
which correspond to ~ 0.45 kpc at the distance of the 
Coma Cluster. 



3.2. Photometric Calibration 

First, we carry out an initial estimation of the pho- 
to metric zero poin t using prev ious photometric catalog 
of IMobasher et all (|2001[ ) and IKomivama et al.l (|2002f) . 
whose observations cover our target fields. However, 
passbands are slightly different between their catalog 
and our data. In order to improve the accuracy of the 
zero point, we re-estimate it using the photometric stars 
from the Sloan Digital Sky Survey (SDSS) DR7 photo- 
metric catalog which covers the Coma Cluster and the 
SDF. This calibr ation method is similar to the one by 
lYagi et all (j2010f ). Photometric zero points of all the im- 
ages are corrected to give the measured counts (ADU) 
consistent with the stellar spectral energy distributions 
(SEDs) given in the Bruzual-Persson-Gunn-Stryker At- 
las 9 which is an exte nsion of the Gunn-Stryker Atlas 
(|Gunn fc Strvkerlll983l ). We select the stars that are de- 
tected in both the SDSS and the Suprime-Cam, and com- 
pare the stellar sequences from our data and the Gunn & 
Strykcr's SED in the two color planes shown in Figure [3l 
Red lines represent the best fits to the Gunn & Stryker's 
stars. Dashed-dotted lines show the best fit to our raw 
data for stars brighter than r = 21, which are calibrated 
based on the initial estimation of zero points. The typi- 
cal number of stars used for the calibration is about 80 
and 30 in the Coma Cluster and the SDF, respectively. 
We correct for these offsets in the zero points by shift- 
ing our stellar sequence to the model sequence computed 
by convolving the Gunn & Stryker's SEDs with the filter 
response curves. The offset value is mostly due to the dif- 
ference of the filter response bet ween the Suprime-Cam 
pas sbands and the fil t ers us ed in IMobasher et al.l (j2001l ) 
and lKomivama et al.l (j2002l ). Specifically, we give the off- 
set values; -0.28, -0.36 and -0.23 mag in f?-band, -0.07, 
-0.08 and -0.11 mag in i?-band of Coma 1, 2 and 3. The 
errors in the final photometric zero points are less than 
0.02 mag and 0.03 mag for the Coma Cluster and the 
SDF images, respectively. 

The Gala ctic extinction i s corr ected using the extinc- 
tion map bv lSchleeel et al.l (|1998l) . Tabled lists the pho- 
tometric zero points and the Galactic extinction in each 
band. 

9 htt p : / / www .stsci.edu/hst/ observatory /cdbs /bpgs . html 



Region Band zp Galactic Extinction 







(mag/ADU) 


(mag) 


Coma 1 


B 


33.04 


0.041 




R 


31.84 


0.025 


Coma 2 


B 


33.85 


0.041 




R 


32.66 


0.025 


Coma 3 


B 


34.23 


0.042 




R 


32.94 


0.026 


SDF 


B 


34.69 


0.078 




R 


34.26 


0.048 


TABLE 2 



List of the photometric zero points and the Galactic 
extinction for the coma cluster and the sdf. 



.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l_ 

r r<22 ~z 
~\ i i i | i i i i | i i i i | i i i i~ 


_ I , , , , I , , , , I , , , 
- ' ».-t= 

r r<21 1 
- , , , , i , , , , i , , , , : 


r<20 " 
: , , , , i , , , , i , , , , i , , , , : 


10 12 3 



-10 12 3 

(fi -r )sDSS 



- 1 1 1 1 1 1 1 - 


.i i i i l i i i i l i i i i l i i i i. 


r r<!22 ~z 
j i i i | i i i i | i i i i | i i i f 


r<21 1 
-, , , , i , , , , i , , , , i , , , , : 




10 12 3 


r<20 " 
: , , , , i , , , , i , , , , i , , , , : 





-10 12 3 

(S -r )sDSS 

Fig. 3. — Color-color plots for stellar objects in the Coma 1 field 
with the initial calibration (black dots). The dot-dashed lines are 
the least square fit to stars brighter than r = 21 and g — r < 1.0 
whose B — g and R — r are calculated from our Suprime-Cam 
measurements and SDSS catalog. The colors of Gunn & Strykcr's 
stars calculated with the filter models are plotted as red points. 
Red lines are the least squares fit to model stars at the range of 
g-r < 1.0. 
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3.3. Object Extraction and Catalog 

We perform object detectio n and photometric mea- 
surement using the SExtractor (|Bertin fc Arnoutslll996l ) 
in double-image mode with the i?-band image as the de- 
tection image. We identify as an object more than five 
connected pixels above 1.5a of the sky background noise. 
We adopt MAG_AUT0 for the total magnitudes and 2" cir- 
cular aperture magnitudes for colors. 

The limiting magnitude i s deriy ed by the method sim- 
ilar to that in lOuchi et al.l (|2004f ) . We measure the sky 
counts in a number of 2" apertures which are distributed 
at random on the image, and obtain the 1 a noise by 
fitting a Gaussian function to the histogram of the sky 
counts. We confirm that there is no failure of a 2" aper- 
ture photometry of detected objects in the i?-band and 
that measurements of i?-band aperture magnitude of the 
object brighter than the limiting magnitude are success- 
ful in double-image mode. The 5a limiting magnitudes 
and the effective areas are summarized in Table [TJ The 
limiting magnitude in each band is found to be ~ 26 mag 
(5a) . Regions around bright stars and giant galaxies (for 
example, NGC 4889, NGC 4874, NGC 4898, NGC 4865, 
NGC 4911 and IC 4051 in the Coma 1), where the de- 
tection of faint objects fails, are masked out. 

The star/galaxy discrimination is performed for ob- 
jects with R < 22. We regard the objects as galax- 
ies which are significantly more extended (FWHM>5.6 
pixels=l."13) than stars (FWHM= 4.9 pixels=0."98) as 
shown in Figure @] As for stellar objects fainter than 
R = 22, the statistical subtraction (for details in ^4.2|) 
is applied in the following analysis. We confirm that the 
LFs remains unchanged if we change the star / galaxy dis- 
crimination limits to R — 21 or R = 23 instead of R = 22. 
Details are given in £15.11 We investigate the contamina- 
tion from globular clusters fainter than R = 22 in detail 
in ^4.41 and it is found to be negligible. 

We evaluate the effective surface brightness of the ob- 
jects (/z e ) in units of mag arcsec -2 . It is defined by the 
mean surface brightness inside the effective radius (ra- 
dius in which half the total flux is emitted) and contains 
some information on the luminosity profile of the objects. 
The distribution of the effective surface brightness as a 
function of total i?-band magnitude is shown in Figure 
[5j The dashed line indicates the sequence of PSF-size 
objects in each field. Saturated objects and stars with 
R < 22 are removed from this figure. 



3.4. Color Classification 

We measure the colors of the objects within 2" aper- 
tures, which corresponds to ~ 0.9 kpc at the Coma Clus- 
ter distance. Figure |5] shows the color-magnitude dia- 
gram (CMD) of the objects in the Coma Cluster and 
SDF. Note that only galaxies arc plotted for R < 22 
where star/galaxy discrimination was applied, while all 
the objects are plotted for R > 22. The color- magnitude 
relation (CMR), indicated by the solid line in each panel, 
is the least-squares fit to all the galaxies in the central 
region that are spectroscopi cally confirmed as mem bers 
of our Coma 1 field from the iMobasher et al.l (|2001l) cat- 
alog. It is represented by 

(B-R) = (-0.029 ± 0M0)R + (1.55 ± 0.21), (1) 



20 
10 

ft 

E 20 
fc 10 




Fig. 4. — FWHM vs. R-magnitudc of all objects (excluding 
saturated objects) in the three observed fields and the SDF. The 
solid line indicates the discrimination between stars and galaxies. 
No star/galaxy discrimination is performed at fainter magnitudes 
(R > 22). 
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Fig. 5. — The effective surface brightness vs. R-magnitude of 
objects in the Coma Cluster and the SDF. The sequence of PSF- 
size objects is shown by the dashed line in each panel. Note that 
we exclude saturated objects and stars brighter than R = 22. 

and broadly consistent with the CMR by (B — R) — 
(-0.045 ± 28)fl + (2.27 ± 0.48) in their bandpass of 
lAdami et al. (20063). 

lAdami et al.l (|2006b| ) reported that most of the faint 
galaxies (21 < R < 25) follow the CMR derived for 
brighter elliptical galaxies (R < 18) in the Coma Clus- 
ter, suggesting that these low-luminosity galaxies ex- 
perienced a passiv e evolution similar to brighter ones. 
lAdami et al.l ((2009) confirmed that the spectroscopically 
identified dwarf galaxies (20 < R < 23) also follow the 
CMR. Assuming that Equation ([TJ is valid down to our 
limiting magnitude (i? ~ 26), we divide the galaxies 
into red and blue populations. We adopt the separation 
boundary which is parallel to the CMR but shifted blue- 
ward by A(B-R) = -0.2, i.e, (B-R) = -0.029E+1.35. 

3.5. Completeness and False Detection Estimate 

Detection completeness of objects in the Suprime-Cam 
images depends on the object magnitudes. We estimate 
the completeness through a series of simulations by de- 
tecting artificial objects in the original images following 
the procedure given in Kashikawa et al. (2004) . Gaus- 
sian profiles with different peaks whose FWHM is equal 
to the seeing size are given to artificial objects. These 
artificial objects are randomly distributed on the object- 
subtracted image, which is -OBJECTS check image made 
by SExtractor. We mask circular areas around objects 
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Fig. 6. — Color magnitude diagram (B — R vs. R) in the Coma 
Cluster and the SDF. The solid line in each panel represents the 
color magnitude relation of bright galaxies in our Coma 1 region. 
We separate red and blue galaxies based on the line shifted blue- 
ward by A(B - R) = -0.2. 

with > 10<7 detection threshold of the sky background 
noise on the image to avoid detections of artificial ob- 
jects embedded above bright objects. The radius of each 
circle is FWHMxl.5 of these bright objects. The artifi- 
cial star magnitudes are distributed every 0.5 mag from 
20.0 mag to 27.0 mag. We extract the artificial objects 
with SExtractor using exactly the same parameters as 
in ^3.3[ and estimate the completeness factor, k(R), i.e., 
the number of the detected objects divided by the num- 
ber of artificial objects in the input list as a function of 
R magnitude. We generate 250 artificial objects for each 
magnitude bin and repeat the process ten times. The 
detection completeness for our catalog is shown by the 
solid curves in Figure [Jj Each image has more than 90% 
completeness at R < 23.5. 

We also show the false detection rate for i?-band de- 
tected objects by the broken curves in Figure[jj It is esti- 
mated by performing photometry on the negative image 
at the same threshold level as the positive i.e. observed 
image using the SExtractor. We confirm that the false 
detection rate is almost 0% down to R = 25 in all fields. 
Therefore, no correction of the false detection rate is per- 
formed in this study. 

We use observed galaxies down to R = 25 [Mr = —10) 
and perform the completeness correction by calculating 
-W (-/^corrected = N (R) ohsQrvod /k(R) for each 0.5 magni- 
tude bin. The magnitude error of artificial objects be- 
tween input and output magnitudes is less than 0.05 mag 
at R = 25 in all images. 

4. CONSTRUCTING LUMINOSITY FUNCTION 

4.1. Number Density 

First, we derive the number density of all the objects 
and red and blue objects as a function of R magnitude 
applying the completeness correction. Figure |8] shows the 
number density of all the objects while Figure [9] shows 
that of red and blue objects in the three fields of the 
Coma Cluster and the SDF. The error bars show the 
Poisson errors. 

4.2. Subtraction of Background/Foreground Galaxies 

The LF is defined as the number density of galaxies per 
unit magnitude interval as a function of absolute magni- 
tude: dN(M) = <j){M)dM 1 where dN(M) is the number 
of galaxies per unit area with magnitudes in the M to 
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Fig. 7. — The detection completeness (solid line) and the false 
detection rate (broken line) for the Coma Cluster and the SDF as 
a function of R-band magnitude. 
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Fig. 8. — The number density of all the objects in the Coma 
Cluster and the SDF plotted against i?-band magnitude. All counts 
are scaled to an area of 1.0 deg 2 and the completeness corrections 
are performed. The Poisson statistical errors are shown as dashed 
lines. 



M + dM range. To derive the intrinsic LFs of cluster 

galaxies, we must correct for the effect of contaminant 

galaxies, which are mainly located behind the cluster. 

We use the statistical background subtraction method 
' n — - — — ~~ 1 ■ — 



e Propris et all 119951 [Bernstein et al.l 11995! 
Phillipps et al.l | 1998t iTrenthaml 119981 : iKambas et al 



20001: lYagi et al.ll2002bl:TPopesso et al.H200a IMilne et alj 



.1. 



We perform the subtraction using the surface 
brightness versus magnitude diagram shown in Figure 
We divide the diagram into meshes with A[i e = 1.0 
mag arcsec~ 2 and Ai? = 0.5 mag and denote the number 
of objects in a mesh as N(/j, e ,R). The number density 
of member galaxies is given by: 

N(jie, R)mb = N(fX e ,R) cl /ki(R) ~ N(pe, R){/k 2 (R), 

where N(p e , R) m h is the number density of member 
galaxies of the cluster per unit area (1.0 deg 2 ), N(fi e , R) c i 
is the number density of all the objects obtained from 
a cluster field, and N(p e ,R){ is the number density of 
all the objects obtained from the SDF, the control field. 
Here, ki(R) and ^(i?) are the completeness factors for 
the cluster field and control field, respectively (Figure [7]). 
This procedure is applied to each of the three Coma fields 
in computing the number density of all the objects and 
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Fig. 9. — The number density of red (top) and blue (bottom) 
objects in the Coma Cluster and the SDF plotted against i?,-band 
magnitude. The color separation is based on the CMR of the Coma 
1 region (for details in All counts are scaled to an area of 1.0 
deg 2 and the completeness corrections are performed. The Poisson 
statistical errors are shown as dashed lines. 

red and blue objects. After the statistical subtraction, 
we summed up N(ji e , i?) m b with respect to \x e to obtain 
the LF of the cluster iV(_R) mb . 

4.3. The Effect of Cosmic Variance 

In order to check the effect of cosmic variance on the 
LF, we re-compute the Coma LF using the SXDS as the 
alternative control field. The SXDS is composed of five 
sub-fields which are referred to as Center, North, South, 
East and West. Each of the SXDS fields has the same 
size of field of view as the SDF and the Coma fields. 
However, the SXDS fields are located at a lower Galactic 
latitude than the SDF and the Coma fields (b ~ 88° for 
Coma, b ~ 83° for SDF, b ~ -60° for SXDS). 

We show in Figure [TU] the LFs in the three regions 
of the Coma Cluster estimated by the statistical back- 
ground subtraction using the the five SXDS fields as the 
control field. The errors of each LF are estimated based 
on Poisson statistics. Figure [TU] indicates that the error 
bars of all LFs overlap in the whole magnitude bins in 
all Coma fields. When the LFs reach the magnitude of 
R ~ —13, their slopes change. The two clear compo- 
nents (bright component and faint component) of differ- 
ent slopes are confirmed in all the LFs. We conclude that 



the effect of cosmic variance is insignificant between the 
SDF and the SXDS fields. In the SXDS-South region, 
some massive structures with masses close to 10 14 Mq 
were discovered by iFinoguenov et ail (|2010D . From their 
results, the SXDS-North region is not polluted by back- 
ground massive structures. We made the LF shown in 
Figure [TU] without masking the massive structure in the 
SXDS-South. We do not see any significant effect of the 
structure. 

In this paper, we adopt the LFs obtained by using the 
SDF as the control field because the SDF and the Coma 
Cluster are located at nearly the same Galactic latitude, 
and the number count of faint stars are similar to the 
two fields. 

4.4. Contribution from Globular Clusters 

At faint magnitudes, we must consider the contami- 
nation from globular clusters (GCs). Unresolved low- 
luminosity galaxies whose angular sizes are similar to the 
seeing size cannot be distinguished from bright GCs that 
appear as point sources at the distance of the Coma Clus- 
ter. One of the funda mental parameters of GC study is 
the specific frequency (jHarris fc van den Berg 



S N = N GC 10' 



0.4(Mv + 15) 



(3) 



where Nqc is the total number of GCs in the galaxy. 
Sm is normaliz ed to a galaxy with an absolute V mag- 
nitude of —15. IMarin-Franch fc Apariciol (|2002f) studied 
the GCs of 17 elliptical galaxies in the Coma Cluster and 
found a mean value of the specific frequency Sm = 5.1. 
On the other hand, the GC luminosity function (GCLF 
hereafter) of a g alaxy is desc ribed by a Gaussian distri- 
bution function (|HarrisHl991[) : 



I \ NgC „-(m-m°) 2 /2a 2 

nQc(m) — — , e v 



(4) 



2na 

where Nqc is the total number of GCs, m° is the 
turnover (peak ) magnitude of the distribution, and a is 
the dispersion. Harris et al.l (|2009D reports the GCLF pa- 
rameters for giant elliptical galaxies in the Coma Cluster 
as V° = 27.7 and a = 1.48. We compute the total GCLF 
around giant red galaxies brighter than Mr — —18 in the 
Coma 1 field based on Equations © and (gj). The GC 
contamination from galaxies fainter than Mr = —18 is 
much smaller than from giant galaxies. No separate esti- 
mation of the total GCLFs is performed in Coma 2 and 
Coma 3 fields since the total GCLFs in these fields are 
supposed to be below that in the Coma 1 field. 



5. RESULTS 



5.1. 



Luminosity Function of All Galaxies 

The top panel of Figure [TTJ shows the LF of all the three 
fields while the bottom panel shows the LFs of Coma 1, 
2, and 3 fields separately. The broken curve in Figure 
[TT1 is the total GCLF of the Coma 1 field computed in 
£14.41 From Figure [TTJ we find that the faintest part of 
the LF (Mr > —12) remains nearly unaffected by GCs 
(< 15%) because the number of faint dwarf galaxies at 
Mr > —12 is about 10 times lar g er tha n that of GCs 
even in the Coma 1. iPeng et "all (|2011| ) observed that 
intracluster GCs in the Coma Cluster have a flat spatial 
distribution within 520 kpc from NGC 4874. The number 
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Fig. 10. — The LFs in the Coma 1 (top-left), the Coma 2 (top-right) and the Coma 3 (bottom) obtained by the statistical background 
subtraction using different control fields, the SDF and the five SXDS fields. Note that each of the five SXDS fields is indicated by center 
(C), north (N), south (S), east (E) and west (W), respectively. The error bars of the LFs come from the Poisson statistics. 

count of intracluster GCs at r — 520 kpc is less than 1% 
of that of GCs around NGC 4874, and it is considered 
to be negligible. Furthermore, we masked the regions 
around giant galaxies on our images as mentioned in ^3-31 
Indeed, at least more than 50% of the total number of 
GCs around these galaxie s are e xcluded from our data, 
if we refer to Harris et "all pOOl . Therefore, the GCLF 
shown by the broken curve in Figure [TT] can be regarded 
as the upper limit of the GC contaminations for both the 
total LF and for each LF of the three fields. 

The top panel shows the clear signature of two compo- 
nents comprising the faint part of the LF. One is domi- 
nant at —19 < Mr < — 13 and the other is dominant at 
— 13 < Mr < — 10. They are hereafter referred to as the 
bright dwarf component (BDC) and the faint dwarf com- 
ponent (FDC), respectively. These distinct components 
appear in all the LFs in spite of the different environ- 
ments such as the Coma 1, 2 and 3. No significant dif- 
ference is found in the shape of the LF among the three 
r egions. 

Mobashcr et al. (2003) reported that the bright part of 
the LF (—23 < Mr < —16) is similar in the central and 
outer regions (includi ng our Coma 2 and C oma 3 fields). 
On the other hand, lAdami et al.l (|2007al ) indicated in 
their Figure 11 that small-scale spatial variations of the 
LF exist in the regions surrounding giant galaxies in the 
center of the Coma Cluster. They divided their observed 
area (~ 0.5 deg 2 ) into 20 subregions of (10' x 10') and 



examined the LFs of the objects in the subregions. In 
this study we are interested in the possible variation of 
the LF over muc h larger spatial scale than that examined 
by lAdami et al.l (|2007al) . i.e., global scale represented in 
terms of the core, intermediate, and outskirt of the clus- 
ter. In this context, we conclude that the shape of the 
LF is very similar in the cluster regardless of the field 
location over the very wide magnitude range so far in- 
vestigated (-23 < M R < -10). 

The ch aracteristic magn itude of the best- fit Schechter 
functio n dSchechterl I1976D t o Coma galaxies is M R ~ 
-20.5 (|Mobasher et alJ l2003). Since most of our galaxies 
are much fainter than this magnitude, we quantify our 
LFs with the logarithmic slope a which is represented, 
as in lAndreon fc Cuillandrd IpOOlL bv 



1 d log ( 



1, 



0.4 8M (5) 

In practice, we derive a by the linear regression method 
as 

log (j> = -0A(a + l)M + const. (6) 

not by fitting the Schechter function to the LF. 

The fitting procedures are performed for the whole 
range of -19 < Mr < -10, and for the BDC and for 
the FDC, separately. We summarize the derived slopes 
a in Table |3l Although slopes of the three fields differ 
slightly, they are consistent within errors. The FDC has a 
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Field 


whole 
-19 < M R < -10 


BDC 
-19 < Mb < -13 


FDC 
-13 < M R < -10 


Coma 1 


-1.77 ± 0.04 


-1.58 ± 0.14 


-3.38 ± 0.28 


Coma 2 


-1.89 ± 0.05 


-1.73 ± 0.22 


-3.58 ± 0.25 


Coma 3 


-1.84 ± 0.07 


-1.54 ± 0.25 


-3.60 ± 0.42 



TABLE 3 

List of the faint-end slope a in each LF for R-band. 
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Fig. 11. — Top: the total J?-band LF of Coma 1, Coma 2 and 
Coma 3. Bottom: the LFs of the Coma 1 (filled circles), Coma 2 
(open circles) and Coma 3 (cross marks) fields, respectively. The 
error bars are based on Poisson statistics. The broken curve indi- 
cates the LF of GCs in the Coma 1 field calcula ted u sing the pre- 
vious results of Marin- Franch & Aparicio (2002) and Harris et al. 
(2009). All counts are scaled to an area of 1.0 deg 2 . Note that 
the LF plots of Coma 2 and Coma 3 are shifted rightward for the 
visibility. 

much steeper slope (a < —3) than the BDC (a ~ —1.6). 

We did not perform star/galaxy separation below the 
limit of R = 22 where many of our galaxies are too small 
to be distinguished from stars or globular clusters (GCs). 
The counts of Galactic stars fainter than this limit in the 
Coma Cluster are subtracted statistically using the cor- 
responding star counts in the SDF. Figure IT21 shows how 
the LF of Coma 1 held changes when we adopt different 
limits. The LFs for R = 21, R = 22 and R = 23 are 
all within the errors. We discuss the contamination from 
GCs for details in 
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5.2. Comparison with LFs from Previous Studies 



Fig. 12. — The ij-band LFs in the Coma 1 (filled circles) when we 
set the star/galaxy discrimination limit to R = 21 (open triangles), 
R = 22 (filled circles) and R = 23 (open crosses). The dashed lines 
indicate errors in Poisson statistics. All counts are scaled to an area 
of 1.0 deg 2 . 

Many observational studies have been performed on 
the Coma Cluster LFs. The LFs of the Coma Cluster in 
the i?-band from previous studies are shown in Figure H31 
Our total LF (red circles) (Coma 1 + Coma 2 + Coma 
3) covers an unprecedentedly wide magnitude range. 

The LF of our BDC agrees with most of the previ- 
ous determinations down to Mr ~ —13, and the LF 
of our FDC is broad l y consistent with the LFs given 
by IMilne et al.l (120071) .iBernstein et all (119951) The LFs 
from lMilne et all (|2007D and iBernstein et al.l (|1995D are 
based on observations of a much smaller region in the 
cluster core. Their effective areas are only < 5% of the 
held of our Coma 1 held. Using our wide and deep 
imaging data, we confirm that the LF in the FDC of 
the Co ma Cluster i s risin g stee ply, similar to t h e hnd - 
ings of IMilne et all (|2007t) and IBernstein et al.l d!995ft . 
at Mr > -13. The plots of lAdami et al.l (|2007bl ) come 
from the data of their held 2 (~ 70 arcmin 2 ) and it is 
inconsistent with our result. 

Their LF have no clear gap which distinguishes be- 
tween the BDC and the FDC, although the small dip is 
found at Mr ~ — 13. However, their LFs in the central 
strip region (~ 300 arcmin 2 ) of the Coma Cluster (in- 
cluding the held 2), seems to be strongly dependent on 
the local environments. Actually, their LF in the held 
2 is very different from that in held 1, next to held 2 
toward west. The LF of their held 1 has deep dip at 
Mr = — 13, and it has no FDC due to the contribution 
from GCs. The difference between their results and ours 
may be due to the difference in the size of observed area, 
i.e., a small area in the cluster core would be suscepti- 
ble to the effect the local environment associated with 
the bright cD galaxy. Our LF is consistent with most of 
the previous results and we show clearly that it is com- 
posed of well-defined two components of the BDC and 
the FDC. 

5.3. Properties of BDC and FDC galaxies 

First, we investigate the color of dwarf galaxies. Figure 
[Ml shows the LF of red and blue galaxies dehned in §3.41 
in the three helds. It is clear from Figure [14] that red 
galaxies are the dominant galaxy population in the BDC 
of all the three helds. The small fraction of blue galaxies 
in the BDC (< 10%) is nearly the same in the three helds. 
On the other hand, the FDC of the Coma 1 is dominated 
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Fig. 13. — The i?-band LFs of galaxies in the Coma Cluster 
from the literature. The red circles indicate our LF in total area 
of Coma 1, Coma 2 and Coma 3. The ordinate shows logarithmic 
counts per 0.5 mag bin per deg 2 . 



by red galaxies, while those of the Coma 2 and 3 contain 
some blue galaxies. The fraction of blue galaxies at the 
faintest magnitude bin (— 11 < Mr < —10) is ~ 10% in 
the Coma 1 field while it is ~ 40% in the Coma 3 field. 

Next, we investigate the relationship between the mag- 
nitude and the surface brightness. Figure [15] shows the 
color-coded density contour map in the magnitude ver- 
sus surface brightness plane. The left, middle, and right 
columns are for the Coma 1, Coma 2 and Coma 3 fields, 
respectively, while the top, middle and bottom lows are 
for all, red and blue galaxies, respectively. The white 
line in Figure [TBI is the line of star-like objects in Figure [5] 
shifted upward by 1.0 mag arcsec -2 , [i e = 0.93Mr + 36.2. 
This line is used here to discriminate between extended 
(low surface brightness) objects and compact (high sur- 
face brightness) objects. We find both extended and 
compact objects in the Coma 1 field, but very few ex- 
tended objects are seen in the Coma 2 and Coma 3 fields. 

Deep spectroscop ic studies of faint galaxies in the 
Com a Cluster (e.g., lAdami et al.1 120091 : iChiboucas et al.1 
l2010T l help us to verify our photometric results de- 
rived from the s tatistical background subtraction. 
lAdami et al.l (|2009l ) presented the faint galaxies (— 14 < 
Mr < —12) which are spectroscopically confirmed clus- 
ter members, follow the CMR of giant member galaxies 
and their surface brightness spans approximately 2 mag 
arcsec -2 at Mr = —13. In addition, they found that the 
compact objects are a minor component which represent 
at most 5% of the member galaxies. This is consistent 
with our results of Coma 1 of Figure [15] The population 
of our compact objects is much smaller than that of ex- 
tended ones down to Mr = — 12. Also IChiboucas et al.1 
(2010) reported the spectroscopic results of dwarf galax- 
ies (—16 < Mr < —11) in the core region of the Coma 1. 
They examined the galaxy membership in terms of their 
properties such as morphology, surface brightness, size, 
and color using 140 dwarf candidates. Their magnitude 
limit is Mr = — 13 for the low surface brightness objects 
and Mr = — 11 for the high surface brightness objects, 
while our limit is Mr = —10 for both of the low and high 
surface brightness objects although we have no spectro- 
scopic confirmation. They obtained the confirmed mem- 
bers of about 50 low surface brightness galaxies (LS- 



BGs) with magnitudes in the range — 16 < Mr < —13, 
and about 20 ultra compact dwarfs (UCDs) in the range 
-13 < Mr < -11. Their results show that the differ- 
ence of the surface brightness between LSBGs and UCDs 
is about 2 mag arcsec -2 at Mr = — 13 (see Figure 13 of 
Chib oucas et al.ll2010T l . We also find that the difference 
of the surface brightness between our extended galaxies 
(objects above the white line in Figure [15]) and compact 
galaxies (objects below the white line in Figure ITB"]) at 
Mr = —13 in the Coma 1 field is about 2 magnitude, 
consistent with their value. Additionally, the magnitude 
distributions of their UCDs are also consistent with our 
compact galaxies. 

Our wide field survey revealed the following. Most 
of BDC galaxies are extended and red. The number of 
such galaxies increases toward the cluster center. Very 
few blue extended galaxies are found in all the three 
fields. On the other hand, FDC galaxies are predomi- 
nantly compact with sizes (FWHM) comparable to the 
seeing size, corresponding to ~ 0.45 kpc at the distance 
of the Coma Cluster. Although most of them are red, 
some are blue. Especially, a lot of blue compact galax- 
ies exist in the FDC of Coma 3 compared to Coma 1. 
In contrast, almost no extended galaxies is found in the 
FDC of Coma 3. The extended galaxies of FDC seem to 
exist only in the dense regions such as the cluster center. 

6. DISCUSSIONS 

First, we compare our total LF of the Coma Cluster 
with those of nearby clusters taken from the literature. 
The slope of the faint part of LF, a, for nearby clusters 
are summarized in Table [4] Note that we re-computed 
the values of a for previous studies by fitting Equation 
@ to the data values given in respective papers, as- 
suming that Hq = 70 km s _1 Mpc -1 , Q, m = 0.3, and 
£1a = 0.7. The slopes were derived by fitting in two 
different magnitude ranges, —15.0 < M < —12.5 and 
-12.5<Af<-10.0 corresponding to the BDC and the 
FDC, respectively. The two fitting ranges are fixed re- 
gardless of the passbands used. We found that the LF 
of the Centaurus Cluster is similar to that of the Coma 
Cluster showing the clear sign of BDC and FDC and 
similar values of the slopes. 

Next, we examine the correlation between the faint- 
end slope (a) and other cluster properties, i.e., the red- 
shift (z), the Bautz-Morgan type, the velocity disper- 
sion (cr), the X-ray temperature (Tx) and luminosity 
(Lx)- The correlation coefficients between a and the 
cluster parameters are listed in Table [5] We evaluate 
the p-value, based on the two-tailed test for the Pear- 
son product-moment correlation coefficient. No signifi- 
cant correlation is confirmed with z, BM type, a, or Lx 
for either of the two components, except for a possibly 
significant correlation with Tx for the FDC (—12.5 < 
M < —10.0). The slope tends to become steeper for 
clusters with higher temperatures. The field LF has a 
much flatter slope (a ~ —1.3) than the slope of the clus- 
ter LFs (e.g..lTrentham et al.l l2005l: iBlanton et all 120051 : 
iLiu et al1l200S)). The X-ray temperature is related to the 
clust er mass plciprich fc Bohringerl 120021 : iStanek et al.l 
[2001) . 

Our comparison and correlation analysis mentioned 
above have led us to the conjecture that the LF in the 
cluster is determined by the intrinsic properties and the 
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Fig. 14. — The red (filled red points) and blue galaxy (open blue points) LFs of i?-band in each field of the Coma Cluster. 



evolutionary stage. Slopes of the BDC and FDC might 
be related to the cluster global mass. In massive clusters, 
there is a tendency for galaxies in the BDC to have a flat- 
ter LF and for galaxies in the FDC to have a steep one. 
Such features of the LFs of high-mass clusters may reflect 
some physical processes that influenced the formation 
and evolution of dwarf galaxies, such as increased intra- 
cluster medium pressure, tidal interaction with crossing 
galaxies, and mergers. 

The cold dark matter scenario assumes that dwarf 
galaxies are formed very e arly in the Univers e . Bas ed on 
cosmological simulatio ns, IRicotti fc Gnedinl (|2005D and 
iBovill fc Ricottil ()2009f ) showed that dwarf galaxies form- 
ing at z > 6 were faint dwarf spheroidal galaxies in the 
present day universe. Their "primordial scenario" has 
dwarf galaxies starting close to their current stellar mass 
of about 10 3 ~ 6 Mq. Some of these galaxies might be 
included in the FDC of the Coma Cluster. 

We found that colors of dwarf galaxies in the Coma 
Cluster seem to be related to the environments. A 
large number of red galaxies are found in the center 
of the Coma Cluster while they get fewer towards out- 
side. The "tidal strippi ng scenario" (|Maver et alj|200lt 
IBournaud fc Dud I2006D predicts that low-luminosity 
galaxies would be formed from more massive gas-rich 
galaxies. In the center of massive clusters, many 
tidally disrupted debris of such gas-rich galaxies are pro- 
duced (|Henriques et al.ll2008l IBournaud et al.ll2008D and 
their gas is effectively removed by ram pressure strip- 



ping fe.g.. iGunn fc Gottl[l972l: iFarouki fc Shapirol[l9& 
[ Takeda et all 11984 iMarcolini et al.l 120031: iGotol 120051 
iMori &: Burkertj|2000f ). This effect may explain the exis- 
tence of many red dwarf galaxies in the center of massive 
clusters. 

We also find a large number of compact galaxies, which 
have comparable sizes to PSF, in all the three observed 
fields. Recent studies suggest the existence of a num- 
ber of UCD galaxies in the core of the C oma Cluster 
(|Chiboucas et al.l l2010t iMadrid et al.|[2010h . It is natu- 
ral to assume that the compact population represented 
in the FDC of the LFs contains the UCDs. Several stud- 
ies suggest that the UCDs are the remnant nuclei of 
dwarf elliptical galaxies that have lost their outer parts 
due to the tidal disruptions during passages close to 
the central clust er galaxy, in a pr o cess called "galaxy 
threshing" (e.g.. jBekki et all l200l iMieske et all 120041 : 
]1P 



IGregg et all [20091. Meanwhile the formation of UCDs 
is linked to the merger of young massive star clusters 
form ed during galaxy merg ers fe.g.. lFellhauer fc Kroupal 
[200llMaraston et al.ll2004D . Thus, major galaxy mergers 
form a large number of compact stellar ob jects that ar e 
progenitors o f GCs and UCDs (|Bournaud et al.l 120081 ). 
IMadrid et all (|2010f ) found that UCD candidates in the 
Coma Cluster are not only red but also blue. If the UCDs 
originate from the interactions of galaxies, the faint end 
of LFs in massive clusters like the Coma Cluster are ex- 
pected to have steep slopes due to the strong tidal field 
and frequent galaxy-galaxy mergers. 
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Fig. 15. — The color-coded number density contours of estimated member galaxies (all galaxies: top row, red galaxies: middle row, 
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in left, middle and right columns, respectively. The statistical background subtraction is carried out here. The number density in each 
mesh is normalized to 1.0 deg 2 and corrected for the detection completeness. The white solid line indicates the discrimination between the 
extended and the compact objects. 



We find that the outskirts contain a large fraction of 
blue galaxie s of —11 < Mr < — 10 compared to the cen- 
tral region. iSmith et al.1 ( 2009f ) found that the passive 
dwarf galaxies are younger in the outskirt than those in 
the core of Coma Cluster, and supported the scenario in 
which many these galaxies are the quenched remnants 
of infalling late-type galaxies. If late-type dwarf galax- 
ies in clusters were star-forming galaxies infalling from 
outside, the steep slopes in the blue galaxy LFs and the 
large number of blue dwarfs in the outskirts of the cluster 
would be explained. 

In summary, the faint part of LF in the Coma Cluster 
seems to consist of some populations of galaxies whose 
origins are different. We suggest that these red galax- 
ies contain primordial dwarf galaxies, UCDs of the tidal 
origin and remnants of disrupted dwarf galaxies, while 
blue dwarf galaxies contain some blue UCDs and infalling 
late- type dwarf galaxies from the outer fields. In order 
to understand their origin, the spectroscopic information 
such as their dynamics and metallicity is essential. Fur- 
ther investigations, we require the spectroscopic observa- 
tions of very faint galaxies down to M = —10 in various 
regions of the Coma Cluster to identify member galaxies. 

7. CONCLUSIONS 

We construct the galaxy luminosity functions (LFs) 
for the Coma Cluster, the richest of the nearby galaxy 
clusters. Deep and wide images of the Coma Cluster 



are obtained with the Suprime-Cam mounted on the 
Subaru Telescope. The resultant LFs cover the range 
of —19 < Mr < —10. Contamination from back- 
ground galaxies is subtracted statistically using the num- 
ber counts of galaxies in a blank field, the Subaru Deep 
Field. 

Our main results are summarized as follows: 

1. We derive LFs in the cluster center region (Coma 1 
field), the NGC 4839 group (Coma 2 field) and the 
outskirts (Coma 3 field). No significant differences 
in the shape of the LFs are found in the different 
fields of the Coma Cluster. 

2. The faint part of LF in the Coma Cluster contains 
two distinct components; the bright dwarf compo- 
nent (BDC; -19 < M R < -13) which has a rela- 
tively flat LF slope and the faint dwarf component 
(FDC; -13 < Mb. < -10) which has a steep slope. 

3. Most of BDC galaxies are extended and red. The 
number of such galaxies increases toward the clus- 
ter center. Very few extended and blue galaxies are 
found in the BDC of all the three fields. The frac- 
tion of blue (both extended and compact) galaxies 
in the BDC (< 10%) is nearly the same in all the 
three fields. 

4. On the other hand, FDC galaxies are predomi- 
nantly compact with sizes comparable to the see- 
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TABLE 4 

List of parameters of nearby clusters. Note that the magnitude is converted into the AB system adopting Ho = 70 km s~ 
Mpc -1 , Q m = 0.3, and f^A = 0.7. The L x and T x from the X-Rays Clusters Database (BAX) are set to Ho = 50 km s _1 

mpc -1 , n m = 1.0. 
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a: velocity dispersion 
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TABLE 5 

List of correlation coefficients between a and the parameter. Note that the p-value under 0.05 is considered significant. 



ing size, FWHM corresponding to ~ 0.45 kpe at 
the distance of the Coma Cluster. Most of them 
are red, but some are blue. Especially, a lot of blue 
compact galaxies exist in the FDC of Coma 3 com- 
pared to Coma 1. In contrast, almost no extended 
galaxies is found in the FDC of Coma 3. The ex- 
tended galaxies of FDC seem to exist only in the 
dense regions such as the cluster center. 
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